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Field of the Invention. 

The present invention relates to a method of determining the digging 

5 position of a digging head in a substantially horizontal digging or excavation and, 
more particularly, to a method that ensures accuracy in determining positions of the 
digging head by lessening the influence of noise magnetic field frequency 
components close to the signal magnetic field to be measured. 
Brief Description of the Prior Art. 

10 A horizontal drilling method, which is one of the boring or trenching 

methods of this kind, uses a small-diameter pipe of 100 mm or less for horizontally 
digging in the ground, and, accordingly, the kind of precision position-determining 
apparatus used in an ordinary small-diameter driving method of excavation cannot 
be placed near a drill. To solve this problem, it is customary in the art to generate 

15 an AC magnetic field by means a coil mounted in the drill head that is detected by 
an above-ground magnetic sensor like a coil to determine the current digging 



This method is simple and easy, but since the magnetic field by the coil is 
a dipolar magnetic field, the field rapidly attenuates with distance from the coil. 
20 The method is defective in that it is unable to achieve a highly reliable 

determination of the digging position when a power line or similar magnetic noise 
source is present in or around the place where the position determination is to be 
performed. 



method that permits highly reliable determination of the digging position by 
detecting the AC signal magnetic field from a coil housed in a drill head even if 
another noise magnetic field is present that might have an effect on the 
determination of the digging position . 



position. 



SUMMARY OF THE INVENTION 



25 



An object of the present invention is to provide a position-determination 



1 



To attain the above object, a method for determining digging position, for a 
boring or trenching excavation, according to the present invention, which (1) 
senses an AC magnetic field from a magnetic field source by means of a magnetic 
sensor provided on the ground and (2) calculates the position of the magnetic field 

5 source from the magnitude and direction of the magnetic field that is sensed. This 
method's construction is characterized by several features. 

In situations where there is, in addition to the signal magnetic field 
generated by the magnetic field source, another noise magnetic field generated by a 
nearby current, at least one of the positions of said magnetic field source - the tilt 

10 angle of said magnetic field source to the vertical direction and the azimuth angle 
of its axial direction of said magnetic field source in a horizontal plane — is 
calculated from a projecting component of the magnetic field sensed by the 
magnetic sensor and projected on a plane or straight line orthogonal to a 
vector-valued direction of the noise magnetic field. 

15 In this regard, it is useful to review the investigation and research relating 

to useful utilization techniques of energy resources entitled "Research for 
Low-Loss Energizing Techniques in Establishment of Advanced 
Telecommunications Networks" (the Composite Development System for New 
Energy Industrial Technique). Here, external noise magnetic fields, which affect 

20 the determination of position in a boring or trenching method of excavation, are 
determined to be primarily generated by a current of some kind. In this case, the 
magnitude of the noise magnetic field varies irregularly with time, but its 
vector-valued direction is constant at each field-sensing position. 

The present invention attains its object through the adoption of the 

25 following steps: 

(A) The direction of the noise magnetic field is detected, and a sensed 
magnetic field in which the noise magnetic field and a signal magnetic field are 
mixed is projected on a plane or straight line orthogonal to the direction of the 
noise magnetic field to obtain a projective component. 



Since the projective component is theoretically free from a component 
derived from the noise magnetic field, at lease one of the position, azimuth angle 
and tilt angle of the magnetic field source is calculated so that the magnitude of the 
projective component (in the case of projection on the straight line) or it magnitude 

5 and direction (in the case of projection on the plane) is substantially equal to a 
theoretically calculated value of a corresponding quantity of a magnetic field 
generated from a magnetic field source, or ,has a minimum difference between the 
former and the latter. The sensed magnetic field is obtained by sensing magnetic 
fields at different positions whose number is determined by how many ones of the 

10 position, azimuth angle and tilt angle of the magnetic field source are unknowns 
and how such unknowns are calculated. 

(B) To obtain the noise magnetic field (1) when the number of noise 
magnetic fields is virtually one, the noise magnetic field is sensed to obtain its 
direction essentially in the absence of a signal magnetic field; or 

15 When the noise magnetic field has frequency components at frequencies 

different from that of the signal magnetic field, the frequency components are 
measured to obtain the direction of the noise magnetic field. 

(2) When the number of magnetic fields is virtually two, the frequency 
components of a first one of the two noise magnetic fields, which are widely 

20 spaced from the frequency components of the second noise magnetic field and the 
signal magnetic field, are measured to obtain the vector direction of the first noise 
magnetic field, and the frequency components of the second noise magnetic field, 
which are widely spaced from the frequency components of the first noise 
magnetic field and the signal magnetic field, are measured to obtain the vector 

25 direction of the second noise magnetic field. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The features of the present invention will be clearly understood from the 
following description taken in conjunction with the accompanying drawings, in 
which: 



Fig. 1 is a perspective view that shows the placement of a magnetic sensor 
according to the method of the present invention; 

Fig. 2 is a vector diagram showing the principle of measurement according 
to the method of the present invention where there is one noise magnetic field 
5 source; 

Fig. 3 is a vector diagram showing the principle of measurement according 
to the method of the present invention where there are two noise magnetic field 
sources; 

Fig. 4 is a flowchart illustrating a procedure for obtaining a projective 
10 component according to the method of the present invention where there is one 
noise magnetic field source; 

Fig. 5 is a flowchart illustrating a procedure for obtaining a projective 
component according to the method of the present invention where there are two 
noise magnetic field sources; 
1 5 Fig. 6 is a flowchart illustrating a procedure for obtaining a signal magnetic 

field where there is one noise magnetic field; 

Fig. 7 is a flowchart illustrating a procedure for obtaining a signal magnetic 
field where there are two noise magnetic fields; 

Fig. 8 is a flowchart showing a measurement procedure according to the 
20 method of the present invention; 

Fig. 9 is a flowchart showing a procedure for calculating the position of a 
signal magnetic field according to the method of the present invention where there 
is one noise magnetic field and where the number of unknowns and the number of 
equations are the same; 
25 Fig. 10 is a flowchart showing a procedure for calculating the position of a 

signal magnetic field source according to the method of the present invention 
where there is one noise magnetic field source and where the number of equations 
is greater than the number of unknowns; 

Fig. 1 1 is a flowchart showing a procedure for calculating the position of a 



signal magnetic field according to the method of the present invention where there 
are two noise magnetic field sources and where the number of unknowns and the 
number of equations are the same; 

Fig. 12 is a flowchart showing the procedure for calculating the position of 
5 a signal magnetic field source according to the method of the present invention 
where there are two noise magnetic field sources and where the number of 
equations is greater than the number of unknowns; 

Fig. 13 is a flowchart showing an example of a procedure for obtaining a 
signal magnetic field by using the magnitude of the signal magnetic field vector; 
10 Fig. 14 is a flowchart showing another example of a procedure for 

obtaining a signal magnetic field by using the magnitude of the signal magnetic 
field vector; 

Fig. 15 is a perspective view showing the placement of a magnetic sensor 
for determining the direction of a noise magnetic field according to the methodof 
15 the present invention when the digging head is the signal magnetic field source and 
is distant from the position of measurement; 

Fig. 16 is a flowchart showing a procedure for calculating the direction of a 
noise magnetic field according to the method of the present invention where the 
digging head is distant from the position of measurement and the signal magnetic 
20 field source stops to generate the magnetic field; 

Fig. 17 is a flowchart showing a procedure for calculating the direction of a 
noise magnetic field according to the method of the present invention when the 
signal magnetic field and the noise magnetic field mix with each other; 

Fig. 18 is a signal frequency spectrum diagram showing how the selection 
25 of a frequency for maximizing a frequency spectrum is used in the process for 
calculating the direction of the noise magnetic field in the flowchart as shown in 
Fig. 13; 

Fig. 19 is a flowchart showing a first method for using the flow of 
candidate vector calculating process in the process of calculating the direction of 



the noise magnetic field in the flowchart as shown in Fig. 17; 

Fig. 20 is a flowchart showing a second method for using the flow of 

candidate vector calculating process in the process of calculating the direction of 

the noise magnetic field in the flowchart as shown in Fig. 17; 
5 Fig. 21 is a flowchart showing a method for using the flow of process of 

evaluating the candidate vector and obtaining the direction of the noise magnetic 

field in the process of calculating the direction of the noise magnetic field in the 

flowchart as shown in Fig. 17; 

Fig. 22 is a signal frequency spectrum diagram showing how the selection 
10 of a frequency for maximizing a frequency spectrum is used in the process of 

calculating the direction of the noise magnetic field in the flowchart as shown in 

Fig. 13; 

Fig. 23 is a signal waveform diagram showing the role of specifying a 
period in which only one noise magnetic field exists by using the fact that the 
15 amplitude of the sensed magnetic field signal diminishes during the OFF period of 
the signal magnetic field that is turned OFF by a predetermined procedure, in the 
process of calculating the direction of the noise magnetic field according to the 
method of the present invention; 

Fig. 24 is a signal waveform diagram showing the role of instantaneous 
20 variations in the amplitude of the sensed magnetic field signal when the signal 
magnetic field is periodically turned OFF in the calculation of the direction of the 
noise magnetic field according to the method of the present invention; 

Fig. 25 is a signal waveform diagram showing the role of instantaneous 
variations in the amplitude of the sensed magnetic field signal when the signal 
25 magnetic field is randomly turned OFF in the calculation of the direction of the 
noise magnetic field according to the method of the present invention; 

Fig. 26 is a flowchart showing a method for calculating the direction of the 
noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, a period 



during which a particular magnetic field is OFF being specified by a sequence that 
starts when a correlation function between the sequence and the sensed magnetic 
field is at its maximum, and the direction of the sensed magnetic field in the 
specified period is regarded as the direction of the noise magnetic field; 
5 Fig. 27 is a flowchart showing a method of calculating the direction of the 

noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, and the 
starting time of a sequence indicating the most likely ON/OFF state of the signal 
magnetic field is calculated from a plurality of times when a correlation function 
10 between the predetermined sequence and the sensed magnetic field is at its 
maximum; 

Fig. 28 is a flowchart showing the method of calculating the direction of 
the noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, where the 

1 5 correlation function between the predetermined sequence and the sensed magnetic 
field is calculated at each of a plurality of time points at which the period of the 
sequence is equally divided and the magnetic field is projected at each point in 
time to a vector formed by the calculated function, and one of vectors whose 
variation is minimal is regarded as the direction of the noise magnetic field; 

20 Fig. 29 is a flowchart showing another method for calculating the direction 

of the noise magnetic field according to the present invention in which the signal 
magnetic field is turned OFF on the basis of a predetermined sequence, where a 
correlation function between the predetermined sequence and the sensed magnetic 
field is calculated at each of a plurality of time points at which the period of the 

25 sequence is equally divided, the magnetic field is projected at each time point to a 
vector formed by the calculated function, and one of vectors whose variation is 
minimal is regarded as the direction of the noise magnetic field; 

Fig. 30 is a flowchart showing a method of calculating the direction of the 
noise magnetic field according to the present invention in which the signal 

7 



magnetic field is turned OFF on the basis of a predetermined sequence, and the 
starting time of a sequence indicating the most likely ON/OFF state of the signal 
magnetic field is calculated from a plurality of points in time when the correlation 
function between the predetermined sequence and the sensed magnetic field is at 

5 its maximum; 

Fig. 31 is a flowchart showing yet another method of calculating the 
direction of the noise magnetic field according to the present invention in which 
the signal magnetic field is turned OFF on the basis of a predetermined sequence, a 
correlation function between the predetermined sequence and the sensed magnetic 

10 field is calculated at each of a plurality of points in time into which the period of 
the sequence is equally divided, the magnetic field is projected at each point in 
time to a vector formed by the calculated function, and one of the vectors whose 
variation is minimal is regarded as the direction of the noise magnetic field; and 
Fig. 32 is a perspective view depicting an example of a magnetic field 

15 sensing frame for use in the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
As depicted in Fig. 1, the position of a digging head 2 under the ground 
surface 1, which creates a signal magnetic field source to be sensed, is determined 
using a magnetic sensor 4 placed on the ground surface 1 at a proper position at the 

20 same time that a power line or similar magnetic noise source 3, which generates a 
noise magnetic field, is placed near the digging position to be determined. More 
specifically, a signal magnetic field vector Hg is provided from the digging head 2 
and a noise magnetic field vector H n is provided from the magnetic noise source 3, 
such as power line. In this case, the magnetic sensor 4 senses a resulting vector 

25 H m that is a combined version of the signal magnetic field vector and the noise 
magnetic field vector H n . 

Referring to Fig. 2, a noise magnetic field of a position vector r and at time 
t is identified as a vector H n (r, t). On the other hand, the signal magnetic field 
generated by magnetic field generating means for position sensing is identified as 

8 



the vector H s (r-r c , 9 C , t). Here the vector 9 C is the angle of orientation of the 
magnetic field generating means, which is defined by three angles of rotation in the 
coordinate system fixed to the ground that is the coordinate system fixed to the 
magnetic field generating means. 

5 Because the noise magnetic field and the signal magnetic field are sensed 

simultaneously, and because the noise magnetic field vector H n ([t] r, t) varies 
randomly with time, it is only possible to extract the signal magnetic field vector 
H s (r-r c , 9 C , t) from the measured magnetic field vector H m (r-r e , 0 C , t) if the 
statistical properties of the noise magnetic field are known and the noise magnetic 

10 field is signal- wise orthogonal to the signal magnetic field. And even if the 
statistical properties separating the noise magnetic field from the signal magnetic 
field are known prior to determining the latter's position, this separation comprises 
a large amount of data, and therefore the conventional scheme is never practical. 
According to the present invention, the direction of a vector e^r) of the 

15 noise magnetic field vector H n (r, t) is obtained by separated means and in a 

coordinate system shown in Fig. 2, where a component vector H m p (r-r c , 9 C , t), is 
projected onto a plane vertical to the direction of a vector e^r) of the measured 
magnetic field vector HJi-i c , 9 C , t), as shown in Fig. 4 (SI, S2, S3). 

20 H m p (r-r c , 0 C , 0=H m (r-r c , 0 C , 0~(H m (r-r c , 0 C , r>e n (r)> n (r) (1) 

This component does not contain the noise magnetic field for the reason given 
below. Since 

25 H m p (r-r c , 0 C , f)=H s (r-r c , 0 C , r)+H n (r, 0 = H s (r-r c , 3, r)+|H n (r, Oh(r). (2) 
it follows that 



H m p (r-r c , 6 C , r)=H s (r-r c , 4, 0"(H s (r-r c , 0 c , 0-e n (r)> n (r). 
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from which it can be seen that the projective component vector H m p (r-r c , 0 C , t) does 
not contain the component of the noise magnetic field vector H n (r, t). 

However, the projective component vector H m p (r-r c , 0 C , t) has lost 

5 information of one axis by the projection onto a plane normal to the direction of the 
vector en(r). That is, since the same projective components are obtained 
regardless of the magnitude of a component parallel to the vector e^r), two 
independent components are obtained. 

Although any method can be used to obtain the two independent 

10 components, it is possible to use the method described below. 

That coordinate axis of a measurement coordinate system C M (described 
below) that is not parallel to the direction e^r) of the noise magnetic field vector 
H n (r, t) is chosen. A unit vector in the direction of the chosen coordinate axis is 
identified as vector em. A vector product, e p> i=e m xe n (r), of the unit vector and the 

15 direction e^r) is perpendicular to the direction of the vector e^r) and hence 
contained in the plane of projection and perpendicular to the coordinate axis em. 
The magnitude of the vector obtained by projecting the projective component 
vector H m P (r-r c , 6 C , t) in the direction of the vector e^i, including the direction of 
the vector, is represented (S4) by a value of Hm/(r-r c , 0 C , t). That is, 

20 

H mfl P (r-r c , 0 c) r)=H/(r-r c , 0 c ,f)e Pil . (4) 
e p ,i = e m xe n (r). (5) 



The next step is to calculate a vector e^ perpendicular to the direction of 
25 the vectors e^i and e^r). The direction of the vector e^ is also perpendicular to 
the direction of the vector e^r) and hence contained in the plane of projection and 
perpendicular to the direction of the vector e^i as well. If the projection of the 
projective component vector H m P (r-r c , 0 C , t) in this direction is represented by 
Hm/(r-r c , 0 C > 0, values Hm/(r--r c , 0 C , t)% yl and Hm/(r-r c , 0 C , t)e^ represent two 

10 



independent vectors into which the projective component vector H m p (r-r c , 9 C , t) is 
separated. Here, 



^ 2 P (r-r c , 0 C , r)=H m p (r-r c , 0 C , f)e M . 
5 e p.2=e p .i x e n (r) 



(6) 
(7) 



By setting the position vector r c and the angle of orientation 0 C of the 
magnetic field source so that a theoretically calculated magnetic field, He(r-r c , 0 C , t) 
generated by the magnetic field source of the position vector r substantially 

10 matches the projective component H/(r-r c , 0 C , t) on the same plane as that of the 
measured magnetic field, it becomes possible to detect the position and orientation 
of the magnetic field source. 

The above analysis treats a situation where the number of noise magnetic 
fields is virtually one. When the number of noise magnetic fields is two, the 

1 5 direction of two noise magnetic fields H^r, t), where j = 1 , 2, may be represented 
by vectors e^O") and e^Cr), and the projection of the measured magnetic field on a 
direction vector e N (r)=eni(r)x e^r) in the coordinate system of Fig. 3 is used as 
depicted in Fig. 5 . That is, 



20 H m p (r-r c , 0 C , 0=H m (r-r c , <?, ,)-(H m (r-r c , <, r)-e N (r)> N (r) (8) 
is calculated (Sll, S12, S13). In the above statement, 



where the symbol "x" indicates a vector product. In this instance, an independent 
component in the projective component is one, that is, only the magnitude of the 
vector. 

If the coordinates of the magnetic field source are represented by vector 



e N (r)=e nl (r)xe n2 (r> 



(9) 



25 



r c (x, y, z) and its angle of orientation by vector 0 C (0 X , 0 y » 9 Z ) where 0 X , 8 y and 9 2 are 
referred to as an angle of rotation, a tilt angle and an azimuth, respectively, when 
the signal magnetic field is axially symmetric, the axis of symmetry is regarded as 
the x-axis, and the angle of orientation is set to vector 0 c (0 y , 0 Z ). 

Further, according to the present invention, it is possible to calculate the 
magnitude H m (r-r c , 0 C ) of each component of the original signal magnetic field 
vector H m (r-r c , 0 C , t) by using the synchronous detection of the measured magnetic 
field vector H m (r-r c , 0 C , t) through the use of a proper one component in Eq. (1) or 
(6) that is a projective component. 

H s (r-r c , 0=(H m (r-r c , 0* m P (r-r c , * e , r)),. (10) 

In the above, H m p (r-r c > 0 C , t) is a proper component of the projective 
component given in Eq. (1) or (6). 

Fig. 4 illustrates a case where the number of noise magnetic field sources is 
essentially one. The processing of Fig. 4 is repeated at each measuring point. In 
a situation where there are essentially two noise magnetic fields, the processing 
shown in Fig. 5 is similarly repeated at each measuring point as is the case with Fig. 
4. 

The magnitude H m (r-r c , 0 C ) of the original signal magnetic field H m (r-r c , 0 C , 
t) can be obtained by calculating Eq. (10) after calculating the projective 
component at each measuring point. Concretely, the process shown in Fig. 6 or 7 
is performed (Sll, S12, S14, S15, S16) or (Sll, S12, S14a, S16). 
HOW TO CALCULATE THE POSITION AND ANGLE OF ORIENTATION 
(1) Definitions of the coordinate system and the angle of orientation: 
A coordinate system is set, one that is fixed to the earth with the z-axis in a 
vertical direction (upward). This system will hereinafter be referred to as the 
measuring coordinate system vector C M . The x- and y-axes are properly set so 
that they form a right-hand system. For example, they are parallel to the direction 
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of projection of the side of a measuring frame is projected onto a horizontal plane. 

The goal here is to calculate the coordinate vector r e (x, y, z) and the 
angle-of-orientation vector 0 C (0 X , 0 y , G z ) of the magnetic field source in the 
coordinate system 

5 Like the coordinate system vector C c of the magnetic field source, a 

coordinate system with the axial direction of the magnetic field source as the x-axis 
is disposed so that the y- and z-axes are horizontal and vertical (upward), 
respectively, when the magnetic field source is placed horizontally. 

The angle-of-orientation vector 0 C of the magnetic field source is defined as 

10 the angle of rotation between the measuring coordinate system vector C M and the 
coordinate system vector C c as described below. First, a coordinate system Ceo 
parallel to the measuring coordinate system vector C M is turned by an azimuth 
angle 0 Z about the z-axis (in either coordinate system). This coordinate system 
will be referred to as coordinate system vector C c i. Next, the coordinate system 

15 vector C ci is turned by a tilt angle 0 y about the y-axis of the coordinate system 
vector Cci itself. This coordinate system will be referred to as the coordinate 
system vector C c2 . Additionally, the coordinate system vector C c2 is turned by an 
angle of rotation 0 Z about the x-axis of the coordinate system vector C c2 itself. 
The angle-of-orientation vector 0 C is determined so that the resulting coordinate 

20 system becomes the coordinate system vector C c . 

(2) Explanation of independent measurands and unknowns: 
When the number of noise magnetic fields is essentially one, the 
measurement of the magnetic field at one place will provide two independent 
measured parameters ("measurands"). When the number of noise magnetic fields 

25 is essentially two, the measurement of the magnetic field at one place will provide 
one independent measurand. In situations where components of the original 
signal magnetic field are obtained by synchronous detection, the measurement of 
the magnetic field at one place will provide three independent measurands. In 
contrast, three coordinate components of the coordinate vector r c (x, y, z) of the 

13 



magnetic field source are unknowns. The azimuth angle 0 2 is not obtainable 
without a different criterion such as the direction of the earth magnetism, and 
therefore the azimuth angle is also an unknown unless the digging head is provided 
with a direction sensor. When the digging head is made of a magnetic material, 

5 or when a buried steel pipe or similar magnetic body lies near the digging head, an 
accurate direction cannot be obtained even if the direction sensor is provided and 
so the azimuth angle 0 2 is an unknown in many cases. The tilt angle 0 y can easily 
be detected by a tilt angle sensor that detects the vertical direction, and hence it is 
known in many cases. The same is true of the angle of rotation 0 X . In particular, 

10 when the signal magnetic field is axially symmetric, if the axis of symmetry is 
regarded as the x-axis, the angle of rotation 0 X becomes meaningless and can be 
ignored. 

At any rate, measured magnetic fields need only to be obtained at different 
positions so that independent measurands are larger in number than unknowns. 
15 Arrangement of the measuring system: 

As shown in Fig. 1, for example, magnetic fields are measured using a 
required number of three-axis magnetic sensors disposed on the ground so that 
their relative positions are known. Because the direction of the noise magnetic 
field changes for each position of measurement, it is necessary to perform for each 
20 position the step (S21) for detecting the position of the noise magnetic field and the 
step (S22) for calculating the projective component of the measured magnetic field 
as shown in Fig. 8 a flowchart of the procedure for determining position according 
to the present invention. Where the component of the original signal magnetic 
field must be calculated by synchronous detection, the processing therefor (in a 
25 specified step) needs to be carried out for each position. 
Flow of measurement processing: 

Detecting (S21) the direction of the noise magnetic field, calculating (S23) 
the position of the signal magnetic field source using the projective component and 
calculating the position f of the signal magnetic field using the signal magnetic 

14 



field component will now be described. 

SITUATION I: Virtually one noise magnetic field: 

When using the projective component when the number of unknowns is 

5 Nu (>1) and the number of noise magnetic fields is virtually one, the projective 
component vector H m p (r-r c , 9 C , t) expressed by Eq. ( 1) is calculated from magnetic 
vectors H m (r-r c , 0 C , t) measured at Nu/2 or more different positions, and the 
position vector r c and the angle-of-orientation vector 0 C are determined (S31) as 
depicted in Fig. 9. In this way, a projective component vector at each magnetic 

10 field sensing point is essentially equal to the projective component vector H/(r-r c , 
0 C , t) of a theoretically calculated magnetic field of the projective component at 
each position of measurement. Thus, it becomes possible to detect (S32) the 
position and orientation of the magnetic field source. 

In a situation where the magnitude of the signal magnetic field 

15 synchronously detected by the projective component is used, the magnetic field is 
measured at Nu/3 or more different positions. When the number Nu(>l) of 
unknowns is even, the number of unknowns and the number of independent 
measurands can be made equal to each other, and if setting magnetic field vector 
H m (r-r c , 0 C , t) is measured at positions N m =Nu/2, a Nu number of such equations as 

20 given below need only to be solved. 

K, q P (r k -r c , 0 C , t)) t -H c , q P (r k -r c> 6 C ) = 0, k = 1, ... ,N m \ q = l,2. (11) 

where q = 1, 2 and represents two directions parallel to the plane of projection and 
25 parallel to each other (S33). Accordingly, H^ q p (r-r c , 0 C , t) and H e>q p (r--r c , 9 C , t), 
where q = 1, 2, represent the magnitude of q-direction components of the measured 
magnetic field vector H m (r-r c , 0 C , t) and the theoretical calculated magnetic field 
vector He(r-r c , 0 C , t), respectively. It is assumed that the vector 0 C represents any 
one of 0 X , 0 y , 0 Z , (0 y , 0 Z ), (0 Z , 0 X ), (9 X , 0 y ), (0 X , 0 y , 0 2 ) and <J>, and that the vector r c 
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represents any one of x, y, z, (y, z), (z, x), (x, y), (x, y, z) and <|), where <() represents 
an empty set. 

For example, when the position vector r c (x, y, z) and azimuth angle 0 Z of 
the magnetic field source are unknown, the magnetic field is to be measured at two 
5 different positions and four equations such as given below are solved, by which it 
is possible to obtain (S34) the position vector r c (x, y, z) and azimuth angle 0 Z of the 
magnetic field source. 

(//^ q P (r k -r c , 0 %% t)) 9 -H %q *(a -r c , 0 2 ) = O, * = 1, 2; q = 1, 2. (12) 

10 

In the above, < . > t represents a time average. 

In Fig. 10 the magnetic field is measured at N m (>Nu/2) different places 
more than Ny/2 in excess of the number Nu of unknowns, since a larger number of 
independent measurands than the number of unknowns can be obtained (S41 , S42, 
15 S43). Here, the position r c (e, y, z) and the angle of orientation 0 Z are calculated 
and yield,: 

c \k =1 J 

where < . > t represents the time average and min { . } means that vectors r c and 8 C 

20 are changed to obtain vectors r c and 0 C , which minimize the contents in { . }. 
Further, a symbol w^q indicates weighting. The Eq. ( 13) can be replaced with the 
following equations (S44). 



H{|lX ^(K q P (r k -r c ^ c) ff) f - H K q p (r k - r c , 6^. (14) 
25 ^njzZ^^/^-r,, 0 C , t]\ -H e , q p (r k -r c , *.f j. (15) 
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It is assumed that vectors r c and 0 C in Eqs. (13), (14), (15) and (16) have the same 

meaning as in the case of Eq. (11). 
5 In this analysis it does not matter whether the measured magnetic field 

vector H m (r-r c , 0 C , t) is a signal that has passed through a band pass filter that 

permits the passage of only those components close to the frequency of the signal 

magnetic field, or a wide-band signal that cannot pass through the band pass filter. 

However, using a signal that has passed through the band pass filter increases the 
10 likelihood that the position of the magnetic field will be determined with a high 

degree of reliability. 

SITUATIQN II: Essentially two noise magnetic fields: 

When the number of unknowns is Nu £1) and the number of noise 

magnetic fields is virtually two, the projective component vector H m p (r-r c , 0 C , t) 
15 expressed by Eq. (4) is calculated from magnetic vectors H m (r-r c , 0 C , t) measured 

at Ny/2 or more different positions, and the position vector r c and the 

angle-of-orientation vector 0 C are determined (S51, S52) as depicted in Fig. 11. 

Thus, the abovementioned projective component vector obtained at each position 

of measurement essentially matches the projective component vector He P (r-r c , 0 C , 
20 t) of a theoretically calculated magnetic field at each position of measurement. 

By this means it is possible to detect (S53) the position and orientation of the 

magnetic field source. 

In this situation, the number Nu (>1) of unknowns and the number of 

independent measurands can always be made to be equal to each other; the 
25 magnetic field vector H m (r-r c , 0 C , t) is measured at different positions of the same 

number as that Nu number of unknowns, and Nu number of such equations given 

below need only to be solved. 
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(tf m P (r k -r c , O c , t)) t -H c ? (r k -r c , 0>O, * = 1, ... 



(17) 



Here, the symbols H^r-rc, 9 C , t) and He P (r-r c , 9 C , t) are the respective magnitudes 
of the projective component vector H m p (r-r c , 9 C , t) and the vector He P (r-r c , 9 C , t) of 
the measured magnetic field vector H m (r-r c , 9 C , t) and the theoretical calculated 
magnetic field vector He(r-r c , 9 C , t). It is assumed that the vector 9 C represents any 
one of 9 X , 9 y , 9 Z , (9 y , 9 Z ), (9 Z , 9 X ), (9 X , 9 y ), (9 X , 9 y , 9 Z ) and <j) and that the vector r c 
represents any one of x, y, z, (y, z), (z, x), (x, y), (x, y, z) and (j>, where § represents 
an empty set. 

For example, when the position vector r c (x, y, z), azimuth angle 9 Z and tilt 
angle 9 y of the magnetic field source are unknowns, the magnetic field is to be 
measured at five different positions and four such equations given below are to be 
solved, by which it is possible to obtain the position vector r c (x, y, z), the azimuth 
angle 9 Z and the tilt angle 9 y of the magnetic field source. 



In the above, < . > t represents a time average, and the vector 9 C = 9 c (9 y , 9 Z ). 

Because a greater number of independent measurands than the number of 
unknowns can be obtained (S61, S62), as seen in Fig. 12, when the magnetic field 
is measured at N m (>Ny) different places more than Nu in excess of the number Nu 
of unknowns, the position r c (e, y, z) and the angle of orientation 9 Z are calculated, 
yielding 



where < . > t represents the time average and min { . } means that vectors r c and 9 C 



K P (r k ^r c , 0 C , r)) f -i/ e p (r k -r c , 0 C )=O, k=h ... ,5. 



(18) 



nun 

r c , a 




(19) 
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are changed to obtain an r c and a 0 C that minimize the contents in { . } . The 
symbol w^q indicates weighting. Eq. (19) can be replaced with the following 
equations (S63). 



^n{Z^|K P ( r ^ r c^ ^ c ,r^-// c P (r k -r c) 0 c fj. 
min|^ w ft ^Jff.'fe - r e . 0 C , ^ -H e p (r k -r c , 0 C )| 



(20) 
(21) 
(22) 



Here it is assumed that r c and 0 C in Eqs. (19), (20), (21) and (22) have the same 

10 meaning as in Eq. (16). 

In the above analysis, it does not matter whether the measured magnetic 
field vector Hn(r-r c , 0 C , t) is a signal that has passed through a band pass filter that 
permits the passage of only those components close to the frequency of the signal 
magnetic field, or a wide-band signal that cannot pass through the band pass filter. 

15 However, the use of the signal that has passed through the band pass filter 

increases the likelihood of determining the position of the magnetic field with a 
great deal of reliability. 

Further, in a situation where the signal magnetic field component is 
obtained by synchronous detection, the direction of the noise magnetic field is 

20 determined and is used to obtain the projective magnetic field. The subsequent 
processing is characteristic of the situations of both one and two noise magnetic 
fields. The projective component vectors H m p (r-r c , 0 C , t) at their respective places 
are calculated from the magnetic field vectors H m (r-r c , 0 C , t) measured at Nu/3 or 
more different places, and a proper one of the calculated vectors is used as a 

25 reference signal to perform the synchronous detection of the measured magnetic 
field vectors H m (r-r c , 0 C , t). Thus, the magnitude H s (r-r c , 0 C ) of the original signal 



magnetic field component is obtained. By determining the position vector r c and 
the angle-of-orientation vector 0 C of the magnetic field source so that the 
magnitude of the original signal magnetic field component and the magnitude 
Hc(r-r c , 6 C ) of the theoretical signal component are equal to each other, it is 
5 possible to detect the position and orientation of the magnetic field source. 

When the number N v (>1) of unknowns is a multiple of 3, the number of 
unknowns and the number of independent measurands can be made to be equal to 
each other; if the magnetic field vector H m (r-r c , 6 C , t) is measured at positions 
N m =Nn/NuV3, Nu number of such equations, as given below, needs only be solved. 

10 

H s (r k -r c , 0 c )-H c (r k -r c , 0 C )=O,* = 1, ... ,N m . (23) 

Fig. 13 shows the flow of this process (S61, S62, S65). Here, Hs(r-r c , 8 C ) and 
Hc(r-r c , 0 C , t) are the signal magnetic field calculated from the measured magnetic 

15 field vectors H m (r-r c , 0 C , t) and the magnitude of the theoretical calculated 
magnetic field vector, respectively. The magnitude Hs(r-r c , 0 C ) of the signal 
magnetic field vector is an averaged quantity, as already explained. It is assumed 
that the vector 0 C represents any one of 0 X , 0 y , 6 Z , (0 y , 6 Z ), (9 2 , 8*), (0„ 0 y ), (0 X , 0 y , 
0 Z ) and <|), and that the vector r c represents any one of x, y, z, (y, z), (z, x), (x, y), (x, 

20 y, z) and <|>, where <j> represents an empty set. 

For example, when the position vector r c (x, y, z) of the magnetic field 
source is an unknown, the magnetic field is to be measured at one place and three 
such equations (as given below) are solved, by which it is possible to obtain the 
position vector r c (x, y, z) and the azimuth angle 0 Z of the magnetic field source. 

25 

H 8 (r-r c )-H c (r-r c )=0. (24) 

The vector r is the vector of the position of measurement. 

Because a larger number of independent measurands than the number of 

20 



unknowns can be obtained when the magnetic field is measured at N m (>Nu/2) 
different places more than Nu in excess of the number N D of unknowns, the 
position r c (e, y, z) and the angle of orientation 9 Z are calculated, yielding 



5 




(25) 



where min { . } means that vectors r c and 9 C are changed to obtain the vectors r c 

and 0 C that minimize the contents in { . }. Further, the symbol w k>q indicates 
weighting. Eq. (25) can be replaced with the following equations. 



Fig. 14 shows the flow of this process (S61, S64, S66). 

DETERMINING THE DIRECTION OF THE NOISE MAGNETIC 

15 FIELD. 

With virtually one noise magnetic field. 

Method (1). The first method for determining the direction of the noise 
magnetic field is used where, in a situation with no signal magnetic field, a noise 
magnetic field is measured using the same measuring system as for measuring the 

20 signal magnetic field. This type of situation is shown in Fig. 15, where the 

detection of the digging position is disturbed by a noise magnetic field source lying 
in the vicinity of the digging route. Another situation is where the signal 
magnetic field source is able to receive a command sent, for example, from the 
ground by some means and is responsive to the command to stop generating the 

25 signal magnetic field. 
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(26) 
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In this instance, as depicted in Fig. 16, if the measured magnetic field is 
represented by a vector H m (r-r c , 6 C , t), because this is essentially a noise magnetic 
field vector Hn(r, t) as shown at a step (S7 1), an average value of its absolute values 
can be used (S72) to calculate the direction vector, en(r)=en, x (r), e^r), e^r)), of 
5 the noise magnetic field as follows: 



«„*(r) = 



K«( r - r c'3,Q|), 

J Z K"( r - r c>*c.4 2 ' 
I a=x, y. z 



a= x, y t z. 



(28) 



Alternatively, a root-mean-square value of the above absolute values can be used to 
10 determine the direction of the noise magnetic field by use of the equation 



J((H mia (r k -r c ,e c , t)j) 

e,a(r)=-r^= a = x,y,z. (29) 

J I ((H„ 8 (r k -r..«.0) ) 

1 a=x, y, z 1 



In this case, the symbol H^r-rc, 0 C , t) is an a component of the measured 
15 magnetic field (the noise magnetic field), and a is any one of x, y and z. 

In the analysis above, it does not matter whether the measured magnetic 
field vector H m (r-r c , 0 C , t) is a signal that has passed through a band pass filter that 
permits the passage of only those components close to the frequency of the signal 
magnetic field, or a wide-band signal that cannot pass through the band pass filter. 
20 However, the use of the signal that has passed through the band pass filter 

increases the likelihood of determining the position of the magnetic field with a 
great deal of reliability. 

Method (2\ 

Step 1 

25 Turning now to Fig. 17, at first the frequency spectrum H m (co) of the 



measured magnetic field vector H m (r-r c , 0 C , t) is calculated (S81, S82, S83) by 
means of the following equation. 

H»=F(H m (r k -r c , r)) (30) 

5 

In this case, the symbol F( . ) represents a Fourier transform; the three components 
x, y and z of the measured magnetic field vector H m (r-r c , 9 C , t) are each Fourier 
transformed. In practice, the above-mentioned frequency spectrum can be 
calculated by FFT (fast Fourier transform) or the like of sampled values of the 
10 measured magnetic field vector H m (r-r c , 9 C , t). 
Step 2 

Next, an angular frequency (ty, where i = 1, 2, . . ., N s , of a large-amplitude 
component, such as a line spectrum, is selected (S84) from the absolute value 
|H m (co)| of the frequency spectrum For the component of each angular frequency 
15 (Oi, where i = 1, 2, . . ., Nns , a candidate unit vector e^r, ©0, where i = 1, 2, . . ., N^, of 
the direction of the noise magnetic field, is calculated (S85) by either method (1) or 
(2) below. 

(1) The absolute values of the Fourier-transformed x, y, and z components 
of the angular frequency concerned are used to calculate the candidate 
20 unit vector e^r, coO, where i = 1, 2, . . ., N^, of the direction of the noise 

magnetic field, by the following procedures: 

(2) 

<H^j |H m >j [H^foj) . 
e nOW = J 1 ; — S i = l, ».,iV ns . (31) 

\a=x,y,z 



25 where H^^Oi), a = x, y, z and i = 1 , 2, . . . , , is a C0i con^onent by the Fourier 
transform of the a conponent of the measured magnetic field. 

(3) A narrow-band filter is formed the pass band of which uses, as the 
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center frequency, the angular frequency CDi, where i = 1, 2, . . ., Nns , and 
the method of Eq. (21) or (22) is used to calculate the candidate unit 
vector en(r, coO, where i = 1 , 2, . . . , . That is, the candidate unit 
vector en(r, 0)0= (e^(r), e„, y (r), ejr)), where i = 1, 2, . . ., ,is 
calculated by 



K & | 7 ^=rff • a= * * * *' = ^ - • (32) 



or by 



, s J((»o,x(r k -«;. 3-^ 0)) 

10 g n , a (r, Q>,) = j 1 , =f=< flf = Jt.y.ci=l.....^- (33) 

J I ((^n,.(r k -r c , ft.q. 0)) 

\<x=x,y t z 1 



Step 3 

The candidate unit vector e„(r, ©0= (e^O"), e^r), e^(r)), where i = 1, 2, . . ., 
^ » is considered as the direction vector e^r) of the noise magnetic field, and for 
1 5 each angular frequency Q>j, where i = 1 , 2, . . . , , the method of Eq. ( 1) is used to 
calculate a projective component vector H m p (r-r c , 0 C , C0i, t). 

H m P (r-r c , <9 C , m,, 0=H m (r-r c , 0 C , t) 

-(H m (r-r c , 0 C , 0-e n (r, ^)> n (r, co), i = l, ...,N m . (34) 

20 

The angular frequency CDj is contained as a variable of the projective component so 
as to explicitly point out that the projective component is dependent on the angular 
frequency 0}, where i = 1 , 2, . . . , N„s . A roper time interval Ttest, which consists of 
Ntest durations Tt, k , where k = 1, 2, . . ., , each having a short time length At, is 
25 chosen, and the variation of the projective component vector H m p (r-r c , 0 C , C0j, t) for 
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each duration T^, where k = 1, 2, . . N^t, is evaluated (S86). It is assumed, here, 
that each duration T^, where k = 1, 2, . . does not overlap other durations. 
More concretely, a variation of statistics v eV ai,k(tf>i) of the N^t, where k = 1, . . ., 
Ntest, which are determined by any one of the methods described below, is 
calculated. 

(1) One or both of the means of absolute values of two orthogonal 
components by 



Veval, M = (K " * • 0 *> ^ 9 = I ^ * = L ~ > *U ' = t - - ^ns" (35) 

where (.) Represents the mean value in the duration and v eva ik is a statistic 
calculated for the duration T^. 

(2) Mean of absolute values 

wM^/fr"* *\s k = X -' N **' i = 1 ' - ■ < 36 ) 

(3) One or both of the means of squares of two orthogonal components by 

v^. lc (<» 1 ) = ((» m . q P (r-r c , 6 e , co„ t)j)^,q = l, 2;k = l, ... , N Kst ; i = \, ... , N m . 

(37) 

(4) One or both of square roots of the means of squares of two orthogonal 
components by 



v eval , >i) = ^((# m . q P (r-r c . 6 e , co„ t)fj^, ? = 1, 2; * = 1, ... , W test ; i = 1 N m . 



(38) 



25 



For the statistics v ev ai,k, where k = 1, . . N^t, calculated by these equations, the 
following equation is applied to obtain (S86) a value of co i)min that is the angular 
frequency 0)j that minimizes var(O)i): 



Jmean k (jv cvalt k (a). )- mean k (v cval> ))Y) 
5 var(^) = J ^ ? ^ / -*J= 1, ... , N m . (39) 

In the above, mean k (.) indicates averaging for the suffix k, that is, 

E() 

mean k (.) = ^ J • (40) 

10 

The magnetic field of the angular frequency co i>min derives from the noise 
magnetic field, and the direction of the noise magnetic field becomes a vector e^r, 

G\min)* 

It should be noted that the angular frequency co i>min , which minimizes 
15 var(CDi), needs only to be measured at one place and need not be obtained at every 
place where to measure the magnetic field. 

With this method, it is also possible to calculate a fluctuation in the 
direction of a vector Hn P (r-r c , 9 C , 0)i, t), as well as an amplitude fluctuation given 
by Eq. (39), and to select the angular frequency co^ at which the direction 
20 fluctuation becomes minimal, or smaller than a predetermined value. In the 
above description, the measured magnetic field vector H m (r-r c , 0 C , t) in Step 1 is a 
wide-band signal, and in Step 3 it does not matter whether the measured magnetic 
field vector H m (r-r c , 0 C , t) is a signal that has passed through a band pass filter that 
permits the passage of only those components close to the frequency of the signal 
25 magnetic field or a wide-band signal that cannot pass through the band pass filter. 
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However, the use of the signal that has passed through the band pass filter 
increases the likelihood of determining the position of the magnetic field with a 
high degree of reliability. 

Fig. 18 shows how to select the frequencies f x (=(i)i/2n) 9 f 2 (=(0 2 /27c), 

5 f n (=(co,/27i, where n=Nns), at which the frequency spectrum becomes maximum; 
Fig. 19 shows the flow of the process, including step S91, for obtaining the 
candidate vector; and Figs. 20 and 21 show the flow of the process, including steps 
S101 and S102 or steps Sill and S112, for evaluating the candidate vector and for 
detecting the direction of the noise magnetic field. 

10 The frequency spectrum H m (a>) need not always be used. That is, the 

signal magnetic field is periodically turned OFF/ON after a predetermined 
procedure; the period Tpenod is divided into equally spaced durations; the candidate 
unit vector e^r, tO, where i = 1 , . . . , Nns, is used in place of the candidate unit vector 
en(r, (£>i)> where i = 1, . . ., N^, which is calculated by Eqs. (25), (26) and (27); and, 

15 thereafter, the duration that minimized the variation by Eq. (33) is calculated by the 
process described above. In this way, the vector e^r, tj) in that duration can be 
adopted as the direction of the noise magnetic field. 

Method (3). In the second method, the large-amplitude angular frequency 
0)i, where i = 1, 2, . . ., N s , is selected from the absolute value |H m (co)| of the 

20 frequency spectrum H m (co). However, the vector e^r, co i>min ) can be obtained as 
the direction of the noise magnetic field in exactly the same manner as in the case 
of the second method, by selecting a proper frequency band near the frequency of 
the signal magnetic field, setting properly spaced test frequencies free from the 
frequency of the signal magnetic field in the frequency band and regarding the test 

25 frequencies as the angular frequency 0)i , as in the second method. 

As is the case with the second method, the angular frequency at 
which var(o)i) becomes minimum needs only to be obtained. 

Fig. 22 shows how to select the frequencies fi(=C0i/27i), f 2 (=C02/27i), . . ., 
fa(=(0)r/27i, where n=Nns) at which the frequency spectrum is maximized. The 
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flow of the subsequent process is the same as depicted in Figs. 17, 19 and 20. 

Method (4). The signal magnetic field is periodically stopped under the 
control of a predetermined procedure. For example, the signal magnetic field is 
periodically stopped at a predetermined time interval. Because the intensity of 

5 the magnetic field being measured decreases while the signal magnetic field is 
stopped, the OFF period of the signal magnetic field is identified by interpreting 
the intensity-decreasing period essentially as a predetermined OFF period. The 
direction of the magnetic field measured during the OFF period is used as the 
direction of the noise magnetic field. The direction of the noise magnetic field 

10 can be obtained using the same method as the first one. Fig. 23 shows how the 
amplitude of the measured magnetic field in this method varies over time. 

Method (5). Here, the signal magnetic field is periodically stopped after a 
predetermined procedure. This is carried out as described below as options (1) 
and (2). 

15 (1) To stop the signal magnetic field on a rectangular-wave-wise: 

The signal magnetic field is repeatedly turned ON and OFF with the period 
Tpenocb for instance. 

=-1, t k +t stop <t<t k+v ' } 

20 

When a sequence s(t) is 1, the signal magnetic field is turned OFF, but when the 
sequence is -1, the signal magnetic field is turned ON. In the above, 



t^-k = T^k=\,2,X... (42) 

25 

(2) To stop the signal magnetic field on a pseudo-random-signal-wise basis, 
for example, when the value is in a random sequence like an M-sequence 
consisting of unit periods of the same length N M , the signal magnetic field is 
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turned ON, but when the value is "1," the signal magnetic field is turned OFF; and 
this sequence is repeated. In this case, the time average of the sequence is set to 
0. 

Fig. 24 shows temporal variations in the amplitude of the measured 
5 magnetic field when the signal magnetic field is stopped by method (1). Fig. 25 
shows temporal variations in the amplitude of the measured magnetic field when 
the signal magnetic field is stopped by method (2). 

Next, as depicted in Fig. 26, the correlation function between the sequence 
s(t) and the norm of the measured magnetic field or the absolute value of its 
1 0 particular component is calculated (S 12 1 , S 1 22, S 1 23) . As the correlation 
function, it is possible to use any one of the equations given below: 

R{T)= j^^lH^r-r,, 0 C , $(t-r)dt. (43) 

R(r) = ^^-^( Hm(r _ rc , 3, ,))%(, _ r y, (44) 

15 R a (T)= £ +Wt ^K«0--i-c. tYt-t)dt, a = x,y,z. (45) 

Kir)= f^VK. OjV*)*. <x = x,y,z. (46) 

In this case, the period for detecting the correlation is set to an integral multiple 
NxTpenod of the period T^od- 

20 The OFF state of the signal magnetic field can be detected (S 1 25) from the 

time x=tsync at which time any one of the above correlation functions is maximized 
(S 124). That is, a sequence s^-t^), which starts at the time T=t sync , is used, and 
when the sequence s(t) is "1," the signal magnetic field is regarded as being OFF; 
in this way, the ON/OFF operation of the signal magnetic field is determined. 

25 In a signal magnetic field OFF period determined in this way, the direction 
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of the measured magnetic field is detected, and that direction is regarded as the 
direction of the noise magnetic field (S126). When the sequence s(t) is (a 
second numerical value), the signal magnetic field is turned ON, and when the 
sequence is "1" (a first numerical value), the signal magnetic field is turned OFF; 
5 this sequence is repeated in this way. 

Method (6). As with the fifth method, the signal magnetic field is 
periodically stopped under control of a predetermined procedure as described 
below. 

(1) To stop the signal magnetic field on a rectangular-wave-wise: 
10 The signal magnetic field is repeatedly turned ON and OFF with the 

period T^od, for example. 



s(t) = h t k <t<t k +t stop . 

(41) 

= -1. '*+'stop^W 



15 When the sequence s(t) is 1, the signal magnetic field is turned OFF. When the 
sequence is -1, the signal magnetic field is turned ON. In the above, 

t M -k = T peAoi ,k = l,%X... (42) 



20 (2) To stop the signal magnetic field on a pseudo-random-signal- wise 

basis: 

For example, when the value is in a random sequence like an 
M-sequence consisting of unit periods T^ of the same length N M , the signal 
magnetic field is turned ON, but when the value is "1," the signal magnetic field is 
25 turned OFF; and this sequence is repeated. 

In this case, the sequence s(t) is chosen so that it changes for each 
predetermined time unit AU. The average of the sequence is "0." 

Next, as depicted in Fig. 27, the correlation function between the sequence 
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s(t) and the norm of the measured magnetic field or the absolute value of its 
particular component is calculated (S 1 3 1 , S 1 32, S 1 33) as in the fifth method, 
possible to use any one of the following correlation functions: 



It is 



R(t) = f;" TW #^ I t)js(t- r)dt. (44) 
f^Ki r - r c 3> * = w- (45) 

faA-^ *)1 <x=*>y>z- (46) 

10 Here, the period for detecting the correlation is set to an integral multiple N T Tperiod 

of the period T^od- 

In this instance, there are present, in general, plural times T=tsync* (k = 1, 2, 

. • • , N^) at which the correlation function is maximized and the maximum value 

exceeds a predetermined value (SI 34). It is assumed that ty^, (k = 1, 2, . . ., 
15 N sync ) is an arrangement of such times in temporal order. When the correlation 

value between the sequence s(t) and the signal magnetic field is appropriate, 

A^sync. k = 'sync, k ~~ 'sync, 1 > & = 2, N sync . (47) 

20 is virtually an integral multiple of the time unit At^. The average of the value 
resulting from the subtraction of an integral multiple Ms^At^ of the time unit 
Atsync, where k = 2, . . ., N synCt k , from At^ , where k = 2, . . ., N^, is then 
calculated as 
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St = J ^ 2 - 

sync 



(^synck-^synck^unit) 



^sync-1 



(48) 



In this case, 

'sync= W+ & S ync ( 49 ) 

5 is the beginning of the sequence signal corresponding to the ON/OFF operation of 
the signal magnetic field. 

Thus, the period during which the signal magnetic field is OFF can easily 
be set based on the sequence s^-t^c). 

By applying the same method as method (1) to the magnetic field vector 
1 0 H m (r-r c , 0 C , t) measured in this period, the direction e^r) of the noise magnetic field 
vector H n (r, t) can be calculated (S138). 

Method 7. With reference to Figs. 28, 29, 30 and, and as is the case with 
the fifth method, the signal magnetic field is periodically stopped, for example, by 
a procedure such as those described below. 
15 (1) To stop the signal magnetic field on a rectangular- wave-wise: 

The signal magnetic field is repeatedly turned ON and OFF with the 
period Tpenod, for instance. 

5(0 = 1, t k <t<t k +t sto? . 

(41) 



20 



25 



and when a sequence s(t) is 1, the signal magnetic field is turned OFF, but when 
the sequence is -1, the signal magnetic field is turned ON. In the above, 

'k + i~'k = W* =1 > 2 ' 3 '"- < 42 ) 

(2) To stop the signal magnetic field on a pseudo-random-signal-wise 
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basis: 

For example, when the value is 1 -1" in a random sequence like an 
M-sequence consisting of unit periods T^t of the same length N M , the signal 
magnetic field is turned ON. When the value is "1," however, the signal 
5 magnetic field is turned OFF; and this sequence is repeated accordingly. In this 
case, the time average of the sequence is "0." 

Next, the correlation function between the sequence s(t) and the measured 
magnetic field H m (r-r c , 9 C , t) is calculated (S141, S142, S143). The period Tpenod 
is divided into equally spaced N^ sections of a length T^, and either the following 
1 0 calculations is performed. 

**(0= { k+rf "km,a(r-r c , <?> *}('-0*.*=i. *=*.y.*- ( 5 °) 

*A) = ('^^(//..aCr-r., O c , t)) 2 s(t-tjdt, k=h ... , A^ v ; a = x,y,z. (51) 



15 Here, the symbol R<x(tk) 9 where a = x, y, z and k = 1, . . N^ , is the time 

correlation between an a component H^r-rc, 0 C , t) of the measured magnetic 
field H m (r-r c , 0 C , t) and the sequence s(t). And 

f k =* 0 +*.r divS k = \... 9 N„ (52) 

20 

The measured magnetic field may also be correlated with the time that is an in- 
multiple of the period Tpenod- That is, 

^(^)={ k+ff,W K a (r-r c , 4,fjs^(f-0**=a *dw; a=x,y,z. (53) 

25 R a (t k > C^^^r-r,, 6>, t)) 2 S mp (t - f k )fc *=1, ... , A^; a= W (54) 
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In this instance, the sequence s(t) is replaced with S^t), where 



^(O^^+^period) 

S mp (t) = s(t), 0<t<r pcriod . 



(55) 



The component B^(T-i Cy 0 C , t k , t) is then calculated, where k = 1 , . . . , , 
of the measured magnetic field H^r-rc, 0 C , t) projected on the vector 



en(tk)=(R x (tk), Ry(tk), Rz(tk)), where k = 1, . . ., NW 



formed by correlation functions R<x(t k ), where oc = x, y, z , corresponding to the 
respective components x, y and z of the measured magnetic field Hm(r-r c , 0 C , t). 
Here, a symbol t k , contained as a variable of the projective component, indicates 
that the projective component depends on the variable t k . 

In this method, the vector e^y can be used as the direction of the noise 
magnetic field through utilization of the time t k . In this time t^ a fluctuation in the 
absolute value of the projective component vector H m (r-r c , G c , t k , t), where k = 1, 
...,Ndi V -- 

Mh> wf^^, ''• k ' 1 - (56) 



- is minimized or becomes smaller than a predetermined value (S 145a). In the 
above, {.) t means the calculation of the time average. 

In addition, variations in the x-component H^ x (r---r c , 9 C , t k , t), y-component 
Hjn, y (r-r c , 0 C , t k , t) and z-component H^r-rc, 9 C > t k , t) of the projective component 
vector H m (r-r c , 0 C , t k , t), where k = 1, . . N^, 
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are calculated, and the vector e^tk) can be used as the direction of the noise 
magnetic field by using the time t k in which the sum of the above-mentioned 
variations 



is minimized or becomes smaller than a predetermined value (S 145b). Here, it is 
possible that the correlation functions R x (tk), R y (tk) and R z 0*) at a certain time t k 
may lose their original signs. Hence, it is necessary to evaluate the fluctuation of 
the projective component at each time t k for four combinations [R x 0*)> RyO*), 
R z (tk)L [Rx(tk), Ry(tk), -R z (tO], [R*(tk), -Ry(tk), Rz(tk)] and [R^, -R y (tk), -R z (tk)]. 

The period during which the signal magnetic field is OFF can also be 
easily set (S145c, S145d) based on the sequence s(t-tk). By applying the same 
method as method (1) to the magnetic field vector H n (r-r c , 0 C , t) as measured in this 
period, the direction e^r ) of the noise magnetic field vector H n (r, t) can be 
obtained (S146). Here, time t k of the equally spaced sections of the time 
length Tdiv derived from the period Tpenod can be used A but it is also possible to use 



Z var A) 



(58) 



or 




(59) 
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the time at which the correlation function given by Eq. (45) or (46) is maximized 
or when the correlation function is maximized and exceeds a predetermined value. 

Determining the direction of the noise magnetic field when the number of 
such fields is virtually two. 
5 Even in a case of two noise magnetic fields, when a first one of the two 

noise magnetic fields has an intensity far greater than that of the second noise 
magnetic field at a first frequency, and the second noise magnetic field has an 
intensity far greater than that of the first noise magnetic field at the second 
frequency, the direction of the first and second noise magnetic field can easily be 
10 calculated through using these frequency components in the measured magnetic 
field. 

If the first or second frequency is close to the frequency of the signal 
magnetic field, the directions of the noise magnetic fields can be determined from 
the measured magnetic field that has passed through a band pass filter that permits 

15 the passage of only those frequencies near the frequencies of the signal magnetic 
field, by the same method as methods (4) or (5), which are used in situations of 
virtually one noise magnetic field. 

When either the first and second frequencies is not equal to the frequency 
of the signal magnetic field, the direction of the respective noise magnetic fields 

20 needs only be calculated by the method used in method (1) for situations of 
virtually one noise magnetic field. 

The present invention can also be put to good use when the signal 
magnetic field generated by the magnetic field source is virtually axially 
symmetric, and the invention permits the determination of position with fewer 

25 magnetic sensors or by a magnetic field being sensed at fewer positions than with 
a magnetic field of low symmetry. 

Further, there are cases in which only one noise magnetic field affects the 
measurement of the digging position and the tilt angle of the magnetic field source, 
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which is an inclination of the axis of symmetry corresponding to the axial direction 
of the signal magnetic field set in the magnetic field set in the magnetic field 
source with respect to the vertical direction, is known. Here, the projective 
component of the magnetic field, measured at each of two or more different 
5 positions and on a plane perpendicular to the direction of the noise magnetic field 
sensed at each of the magnetic field sensing positions, is calculated. The position of 
the magnetic field source and its azimuth angle (the direction of the axis of 
symmetry in a horizontal plane) can be calculated from the above projective 
component. 

10 In addition, according to this method of the subject invention, when 

virtually one noise magnetic field affects the measurement of the digging position, 
one can calculate the projective component of the magnetic field, measured at each 
of three or more different positions, on a plane perpendicular to the direction of the 
noise magnetic field sensed at each of the magnetic field sensing positions. The 

15 position of the magnetic field source, its tilt angle (an inclination of the axis of 
symmetry corresponding to the axial direction of the signal magnetic field set in 
the magnetic field source with respect to the vertical direction), and the azimuth 
angle of the magnetic field source (the direction of the axis of symmetry in the 
horizontal plane) can be calculated from the above projective component. 

20 There are situations in which virtually two noise magnetic fields alone 

affect the measurement of the digging position and the tilt angle of the magnetic 
field source, which is an inclination of the axis of symmetry corresponding to the 
axial direction of the signal magnetic field set in the magnetic field source with 
respect to the vertical direction, is known. Here, the projective component of the 

25 magnetic field, measured at four or more positions, on a straight line perpendicular 
to both the direction of a first of two noise magnetic fields sensed at each magnetic 
field sensing position and the direction of the other of the two noise magnetic 
fields sensed at the same position, can be calculated, and the position of the 
magnetic field source and its azimuth angle that is the direction of the axis of 
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symmetry in the horizontal plane can be calculated from the above projective 
component. 

Further, in situations where virtually two noise magnetic fields alone affect 
the measurement of the digging position, the projective component of the magnetic 

5 field, measured at five or more positions, on a straight line perpendicular to both 
the direction of the first of two noise magnetic fields sensed at each magnetic field 
sensing position and the direction of the other noise magnetic field sensed at the 
same position can be calculated. Here, the position of the magnetic field source, 
its tilt angle that is an inclination of the axis of symmetry corresponding to the axial 

10 direction of the signal magnetic field set in the magnetic field source with respect 
to the vertical direction, and the azimuth angle that is the direction of the axis of 
symmetry in the horizontal plane can be calculated from the above projective 
component. 

Moreover, when virtually two noise magnetic fields alone affect the 
15 measurement of the digging position, the frequency component of the first of two 
noise magnetic fields, near which the remaining second noise magnetic field and 
the signal magnetic field have substantially no frequency components, is measured 
in order to permit detection of the direction of the first noise magnetic field in 
terms of a vector. Here, as well, the frequency component of the second noise 
20 magnetic fields, near which the first noise magnetic field and the signal magnetic 
field have substantially no frequency components, is measured in order to detect 
the direction of the second noise magnetic field in terms of a vector. 

In the subject invention, a three-axis magnetic sensor that senses three 
magnetic fields orthogonal to one another at substantially the same position is 
25 effective as a magnetic sensor. 

However, the magnetic sensor in the present invention may be any kind 
sensor as long as it can function in three magnetic fields orthogonal to one another 
at substantially the same position. The three-axis magnetic sensor is suitable 
because it senses three magnetic fields orthogonal to one another at substantially 
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the same position. In alternatively, one magnetic sensor that is capable of sensing 
a magnetic field in only one direction can be turned in the same position to face 
three different orthogonal directions in sequence to sense the three orthogonal 
magnetic fields. 

5 In employing the present invention, one can use a frame such as frame 12 

depicted in Fig. 32, which has a magnetic sensor fixing means 1 1 mounted on it so 
that it can fix three-axis magnetic sensors and a tilt angle gauge 1 3 for detecting the 
inclination of the frame with respect to a vertical direction. The position of each 
magnetic sensor fixing means on the frame is known, and the magnetic sensor 

10 fixing means acts to fix the magnetic sensor in a predetermined orientation to the 
frame. The magnetic sensor fixing means 1 1 can be provided with three faces 
orthogonal to one another, and has a mechanism that fixes the magnetic sensor at a 
predetermined angle when a predetermined face of the sensor case is pressed 
against any of the three faces of the sensor fixing means. One magnetic sensor or 

15 one three-axis magnetic sensor is fixed to these magnetic sensor fixing means one 
after another to sense the magnetic fields. 

In another alternative, a plurality of magnetic sensors may be fixed in 
predetermined orientations to the frame 12 at a plurality of positions to 
simultaneously sense magnetic fields at this plurality of positions. 

20 As described above, the present invention employs a frame provided with a 

plurality of magnetic sensor fixing means each capable of removably or fixedly 
mounting a three-axis magnetic sensor and a tilt angle sensor capable of detecting 
the tilt angle of an orthogonal coordinate system to the vertical direction. The 
magnetic fixing means are mounted on the frame so that their positions and 

25 orientations are known. And the magnetic sensor is removably or fixedly 

mounted on a required number of magnetic sensor fixing means to sense magnetic 
fields, while the tilt angle of the frame during magnetic field sensing and the 
orientation of the magnetic sensor at each magnetic sensor mounting position 
toward the frame are used to calculate, from the magnetic field sensed at each 
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magnetic sensor mounting position, the sensed magnetic field, a noise magnetic 
field and a signal magnetic field as vectors in a coordinate system fixed to the 
ground. 

The magnetic field generating means for the signal magnetic field in the 
5 present invention may be a coil. The magnetic field generating means may also 
be one electric wire, or one electric wire that is straight only near the place of the 
determination of position. 

Even if a buried power line, railroad tracks, or similar noise magnetic 
sources are present near a construction site, the subject invention will afford a 
10 highly reliable determination of position and will note be affected by the noise 
magnetic fields generated by such noise magnetic field sources. 

The present invention is intended for measuring the digging position in the 
non-open-cut method of excavation, but is applicable as well to many technical 
fields that determine position by means of sensing magnetic fields. 
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